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Abstract 

We perform a comprehensive analysis of the Minimal Supersymmetric Standard 
Model (MSSM) in the scenario where the scalar partners of the fermions and the Higgs 
particles (except for the Standard-Model-like one) are assumed to be very heavy and 
are removed from the low-energy spectrum. We first summarize our determination of 
the mass spectrum, in which we include the one-loop radiative corrections and resum to 
all orders the leading logarithms of the large scalar masses, and describe the implemen- 
tation of these features in the Fortran code SuSpect which calculates the masses and 
couplings of the MSSM particles. We then study in detail the phenomenology of the 
model in scenarios where the gaugino mass parameters are non-universal at the GUT 
scale, which leads to very interesting features that are not present in the widely studied 
case of universal gaugino mass parameters. We discuss the constraints from collider 
searches and high-precision measurements, the cosmological constraints on the relic 
abundance of the neutralino candidate for the Dark Matter in the Universe - where 
new and interesting channels for neutralino annihilation appear - and the gluino life- 
time. We then analyze, in the case of non-universal gaugino masses, the decays of 
the Higgs boson (in particular decays into and contributions of SUSY particles), of 
charginos and neutralinos (in particular decays into Higgs bosons and photons) and of 
gluinos, and highlight the differences from the case of universal gaugino masses. 



1 Introduction 



The main reason for introducing low energy Supersymmetry (SUSY) [1] in particle physics 
was its ability to solve the naturalness and hierarchy problems [2]. Indeed, SUSY prevents the 
Higgs boson mass from acquiring very large radiative corrections: the quadratically divergent 
loop contributions of the Standard Model (SM) particles to the squared Higgs masses are 
exactly canceled by the corresponding loop contributions of their supersymmetric partners. 
This cancellation stabilizes the huge hierarchy between the Grand Unification (GUT) and the 
electroweak symmetry breaking (EWSB) scales and no extreme fine tuning of parameters 
relevant to the mechanism of EWSB is required for this purpose. However, SUSY is not 
an exact symmetry and there is still a residual contribution to the Higgs masses that is 
proportional to the mass differences between the SM particles and the new SUSY particles. 
Therefore, for the fine-tuning problem not to be reintroduced in the theory, the mass Ms 
of the new particles should be at most of the order of the TeV. The requirement of no 
fine tuning is the main reason for expecting SUSY particles to be accessible at the next 
generation of high-energy colliders, in particular at the Large Hadron Collider (LHC) and 
the International Linear Collider (ILC). 

Nevertheless, there is no compelling criterion to define the maximal acceptable amount 
of fine tuning [3] and the choice of the upper bound on the SUSY scale Ms is somewhat 
subjective. Thus, it might well be that Ms is significantly larger than one TeV, in which 
case the SUSY particles would be too heavy and not observable at the LHC. However, there 
are two more important motivations for SUSY that also call for some light SUSY particles: 
the consistent unification of the three gauge coupling constants at the GUT scale [4] and the 
solution to the Dark Matter (DM) problem [5]. Indeed, the SM slope of the evolution of the 
three gauge couplings has to be modified early enough by some SUSY particle contributions 
to achieve unification, and the DM problem calls for the existence of a new stable, neutral and 
weakly interacting particle that is not too heavy in order to have the required cosmological 
relic density. However, it has been pointed out [6-8] that, for SUSY to provide solutions to 
the unification and DM problems, only gauginos and higgsinos, the spin-| superpartners of 
the gauge and Higgs bosons, need to be relatively light, with masses of the order of the EWSB 
scale. The scalar partners of the SM fermions sit in complete irreducible representations of 
SU(5), therefore they could be very heavy without spoiling gauge coupling unification. For 
Ms ^ 1 TeV the model would be extremely fine-tuned, though, and one would have to give 
up SUSY as the solution to the hierarchy problem. 

An interesting feature of such a scenario, commonly known as Split Supersymmetry 
[6,7], is that it is much more predictive than the general Minimal Supersymmetric Standard 
Model (MSSM) [9]. Indeed, it is well known that, in the most general case, the MSSM 
has a very large number of free parameters, (9(100), most of them related to the sfermion 
sector. Even if one constrains the model to have a viable phenomenology, as in the so- 
called phenomenological MSSM [10] where minimal fiavor mixing and CP violation as well 
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as universality of the first- and second-generation sfermion masses are assumed, there are still 
more than 20 free parameters left. On the other hand, if the masses of all the scalars (except 
for one SM-like Higgs doublet) are pushed to large values only a handful of parameters 
are needed to describe the low-energy SUSY theory. As a by-product, such an assumption 
cures many problems of the general MSSM (such as the excess of flavor and CP violation, 
fast proton decay with dimension-5 operators, etc.) while still allowing for gauge coupling 
unification and a good DM candidate, the lightest of the neutralinos. 

Indeed, in the MSSM with heavy scalars, besides a common value Ms of the soft SUSY- 
breaking sfermion mass parameters, the basic inputs are essentially the three gaugino masses 
Mi,M2,M3, which can be unified to a common value at the high-energy GUT scale, the 
Higgs-higgsino mass parameter /x which is no longer fixed by the requirement of proper 
EWSB as in the general MSSM, and the parameter tan/5. The trilinear Higgs-sfermion 
couplings are forced to be small by the same symmetry that protects the higgsino and 
gaugino masses, and they play a very minor role. One can then derive in an exhaustive 
manner the relationship between this small number of inputs and the physical parameters, 
i.e. the (super)particle masses and couplings. However, if the scalars are very heavy, they will 
lead to significant quantum corrections in the Higgs and gaugino-higgsino sectors, enhanced 
by potentially large logarithms of the ratio between the EWSB scale and the scalar mass 
scale, log(MEwsB/^5)- In order to have reliable predictions, one has to properly decouple 
the heavy states from the low-energy theory and resum the large logarithmic corrections by 
means of Renormalization Group Equations (RGEs). 

In fact, from a more practical point of view, in most (if not all) of the numerical RGE 
codes that calculate the Higgs and SUSY particle spectrum of the MSSM [11,12], one cannot 
assume too large masses for the sfermions, Ms > a few TeV, as the programs become 
unreliable. Indeed, one has to perform major modifications to the general algorithms of such 
programs, which involve three main steps: the RG evolution, the consistent implementation 
of EWSB and the calculation of the mass spectra including the radiative corrections. In 
particular, one needs to properly decouple the heavy scalars and include the intermediate 
scale Ms in the RG running of parameters back and forth between the low-energy scales Mz 
and Mewsb and the high-energy scale Mqut! ignore the implementation of (radiative) EWSB 
as a large fine tuning is tolerated; resum the large logarithms in the radiative corrections to 
the calculated physical masses and couplings as mentioned above. 

The purpose of this paper is to analyze, in the most comprehensive way, the SUSY and 
Higgs particle spectrum in the MSSM with heavy scalars, with Ms ranging from a few TeV 
to ~ 10^^ GeV. We start by summarizing our determination of the particle spectrum of the 
low-energy theory in which the heavy scalars have been integrated out. We determine the 
masses and couplings of the Higgs boson and of the higgsinos and gauginos at the weak scale 
by solving the one-loop RGEs of Spht SUSY [7, 13], thereby resumming to all orders the 
leading corrections that involve logarithms of the large scalar masses. We also include one- 
loop radiative corrections to the Higgs boson, chargino and neutralino masses, by adapting 
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to the Split-SUSY case the MSSM formulae presented in ref. [14]; we investigate the variation 
of the masses with the renormahzation scale, which can be considered as a rough estimate 
of the higher-order corrections. We also describe the implementation of these features in 
SuSpect [11], one of the public computer codes that calculate the MSSM mass spectrum. 

We then study the regions in the parameter space that are compatible with the presently 
available collider data as well as with the constraints on the relic density of DM and on 
the gluino lifetime. In the former case we reemphasize that a new possibility [15], which is 
not present in the usual MSSM, has to be considered: the efficient annihilation of the DM 
lightest neutralinos through the exchange of the Higgs boson which decays into a real and 
a virtual W boson, the latter subsequently decaying into two massless fermions. Finally, we 
have adapted to the Split-SUSY scenario the programs HDECAY [16] and SDECAY [17] for the 
calculation of the decay widths and branching ratios of the MSSM Higgs bosons and of the 
SUSY particles, respectively. We briefly discuss some features in the decays of the Higgs 
boson and the gauginos, focusing on interesting channels such as Higgs decays into invisible 
neutralinos and into two photons, chargino/neutralino decays into lighter ones and a Higgs 
boson and a photon, and gluino decays into a gluon and a neutralino. 

Since the number of input parameters in the low-energy theory is rather limited, one can 
relax the usual assumption [18] of unified gaugino masses at the GUT scale. Mi = M2 = 
M3 = mi/2, and still have a rather predictive model. Therefore, for the sake of generality, we 
will also consider specific patterns of non-universal gaugino masses [19] and, for illustration, 
we will discuss two different models: one in which SUSY-breaking occurs via an F-term 
that is not a SU(5) singlet [20] and another based on an orbifold string model [21]. In these 
models, the gaugino mass parameters at the electroweak scale can be widely different from 
the pattern of the universal scenario. In particular, the neutralino that is the lightest SUSY 
particle (LSP) can be, most of the time, either wino-like (as in anomaly mediated SUSY- 
breaking models for instance) or higgsino-like (which also occurs in the universal scenario 
for small values of the parameter fi), leading to a near degeneracy of the LSP neutralino 
mass with the lightest chargino mass. In addition, in some of these scenarios, the LSP can 
also be close in mass with the gluino and/or a bino with a very small mass. This leads to 
several interesting phenomenological features that we will discuss in some detail, such as 
rapid annihilation of neutralino DM through the exchange of a Z boson or co-annihilation 
of the gluinos which can be the next-to-lightest SUSY particle (NLSP). 

The rest of the paper is organized as follows. In the next section we summarize the model 
with heavy scalars and discuss the implementation of the RGEs, the radiative corrections to 
the Higgs mass as well as to the gluino, chargino and neutralino masses, and the boundary 
conditions on the soft SUSY-breaking gaugino masses. In section 3 we summarize the various 
constraints on the model parameter space from collider data and cosmology. In section 4 we 
present some results on the Higgs boson decays involving SUSY particles and the decays of 
the chargino and neutralino states as well as the gluinos. Conclusions are given in section 5. 
Finally, a set of useful formulae is collected in the appendix. 
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2 The low-energy effective theory 



In this section we summarize our determination of the mass spectrum of the effective theory 
vahd below the scale Ms at which all the heavy scalars of the MSSM, i.e. the sfermions and 
one Higgs doublet, are removed from the spectrum. We also discuss the different boundary 
conditions on the soft SUSY-breaking gaugino mass parameters that will be considered in 
the phenomeno logical analysis of sections 3 and 4. Finally, a brief summary is given of how 
the model is implemented in the Fortran code SuSpect. 



2.1 Determination of the mass spectrum 



If the common mass of the scalars is assumed to be very large, Ms ^ 1 TeV, the low-energy 
theory contains, besides the SM particles, the higgsinos H^, the gluino g, the wino W 
and the bino B. Omitting the gauge-invariant kinetic terms, as well as non-renormalizable 
operators suppressed by powers of the heavy scale M5, the Lagrangian of the effective theory 
reads [7] 

C D m^H^H - ^ {H^Hf - \^hlqjU,eH* + hffljdiH + hl/jdH 
+ ^M3^V + WW + ^Mi BB + fiH^eHd 



+ ^[g^a''W''+~giB)H^ 



-~gda''W'' + g:,B] Ha + h.c. 



where a"" are the Pauli matrices, e = ia^ and i,j are generation indices. The SM-like Higgs 
doublet H is a linear combination of the two MSSM Higgs doublets Hu and Ha, fine-tuned 
to have a small mass term m^: 



H = — cos P e H*, + sin j3 Hu 



(2) 



At the high scale Ms the boundary conditions on the quartic Higgs coupling and on the 
Higgs-higgsino-gaugino couplings of the effective theory are determined by Supersymmetry: 

1 



\{Ms) 
9u{Ms) 



[g^Ms) + g'^Ms)] cos^ 2/3 + A^hA 



g (Ms) sin p , 
g'{Ms) sin(3 



UMs)=g (Ms) cos (3 , 
g^iiMs) =g' (Ms) cos (3 . 



(3) 

(4) 
(5) 



where g and g' are the SU(2) and U(l) gauge couplings. Note that tan P is not a parameter of 
the low-energy effective theory, and it enters only the boundary conditions on the couplings. 
In fact, contrary to what happens in the usual MSSM, tan/3 is not interpreted here as the 
ratio of two Higgs vacuum expectation values, but rather as the fine-tuned angle that rotates 
the two Higgs doublets into one heavy and one light, SM-like doublet. 
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In the boundary condition to the quartic Higgs couphng, we include also a one-loop 
threshold correction of 0{hf), where ht = rrit/v is the top quark Yukawa coupling, 



Xt = At — fi/ tan/3 is the left-right mixing of the stop squarks, with At being the trilinear 
Higgs-stop coupling; here, Ms is interpreted as the average of the two stop masses. Note, 
however, that in Split-SUSY models the trilinear coupling At cannot be too large, since it is 
protected by the same symmetry that keeps the gaugino and the higgsino mass parameters 
small. The threshold correction in eq. ([6]) can thus be relevant only for relatively small values 
of Ms- Beyond tree level the boundary conditions in eqs. ©-(IS]) are valid only in the DR 
renormalization scheme. The one-loop electroweak corrections to eqs. (I3])-(I5]) that account 
for the shift to the MS scheme are given in the appendix. 

The gauge and third-family Yukawa couplings as well as the vacuum expectation value 
(vev) of the SM Higgs field (normalized as u ~ 174 GeV) are extracted from the following set 
of physical inputs [22,23]: the strong gauge coupling as{Mz) = 0.1187; the electromagnetic 
coupling a{Mz) = 1/127.918; the Z-boson mass Mz = 91.1876 GeV; the Fermi constant 
Gf = 1.16637 X 10"^ GeV-2; the physical top and tau masses Mt = 170.9 ± 1.8 GeV and 
Mr = 1.777 GeV and the running bottom mass rrihimi,) = 4.25 GeV. We use one-loop 
formulae from ref. [14] to convert all the physical inputs into running parameters evaluated 
in the MS scheme at a reference scale that we choose as equal to Mz- To this purpose we 
need to adapt to the Split-SUSY scenario the formulae of ref. [14], which were originally 
derived for the MSSM in the DR scheme: we remove the contributions of the heavy scalars 
and insert appropriate DR-MS conversion factors. A summary of the relevant formulae is 
given in the appendix. 

The gaugino masses are given as input at the GUT scale, defined as the scale where 
the two gauge couplings gi = a/5/3 and g2 = g unify, and evolved down to the scale 
Ms by means of the one-loop RGEs of the MSSM. In addition to the minimal case where 
the three gaugino masses unify at the GUT scale, i.e. Mj(MGUT) = ""^1/2; our analysis will 
consider scenarios of SUSY breaking in which the boundary conditions at the GUT scale are 
different. Finally, the n parameter has to be provided as an independent input, contrary to 
the constrained MSSM case in which it can be extracted from the EWSB conditions, and 
we take it as a running parameter evaluated at the scale Mz- 

The parameters of the Lagrangian in eq. ([1]) are then evolved to a common renormal- 
ization scale Qw of the order of the weak scale, by means of the one-loop RGEs of the 
Split-SUSY model, which we take from ref. [7]. Since some of the boundary conditions on 
the parameters are given at the SUSY scale Ms and the others are given at the weak scale 
Mz an iterative procedure is necessary. The resulting couplings evaluated at the weak scale 
Qw account for the all-order resummation of the leading logarithmic corrections involving 
powers of \og{Ms/Qw)- 




5 



Once the iteration for the determination of the Lagrangian parameters converges, the 
physical masses of the Higgs boson, the charginos and the neutrahnos are computed at the 
scale Qw including one-loop radiative corrections; the gluino mass is computed separately 
at the scale M3. The relation between the physical Higgs boson mass and the quartic 
coupling A computed at the generic scale Q reads 



Mh = \I^^ [1 + S'''{Q) + SHQ)] ■ (7) 
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The SM contribution can be found in ref. [24], and contains terms enhanced by M'^ 
coming from top-quark loops. The remaining term 6^ is the contribution of the diagrams 
involving charginos and neutralinos, and reads 



5^ = - 
2 



nl^iml) ^ n^^(o) 



where m\j = 2 Xv"^ and = ^g^v^ are tree-level Higgs and W masses, while T^, H^jj and 
^ww denote the chargino and neutralino contributions to the Higgs boson tadpole, the Higgs 
boson self-energy and the W boson self-energy at zero external momentum, respectively. The 
explicit dependence of tadpole and self-energies on the renormalization scale compensates, 
up to higher-order (i.e. two-loop) effects, the implicit scale dependence of the Higgs quartic 
coupling in eq. ([7]). We extracted the formulae for the Higgs tadpole and the self-energies 
from the MSSM results of ref. [14], by appropriately rotating the Higgs fields, dropping the 
contributions of heavy scalars and expressing the couplings of charginos and neutralinos in 
terms of the effective Higgs-higgsino-gaugino couplings gu,d and We have also checked 
that our results for 6^ agree with the recent computation of ref. [25] . The relevant formulae 
can be found in the appendix. 

The logarithmic dependence of Mh on the sfermion mass scale Ms has been discussed 
in earlier papers [6,7,13,25], but we show it for completeness in fig. [1], choosing the SUSY 
input parameters as At = 0, fi{Mz) = rni/2 = 500 GeV and tan/3 = 2 or 30. The top 
Yukawa coupling entering the one- loop corrections to the Higgs mass, eqs. (EI) and (JTj), is 
expressed in terms of the running top mass, which is in turn extracted from the physical top 
mass via suitable threshold corrections. For each choice of tan /3 the solid lines in fig. [1] are 
obtained with the central value of the physical top mass, Mt = 170.9 GeV, and the dashed 
lines correspond to a ±lcr variation in the top mass. It can be seen that, for relatively low 
values of Ms, going from small to large tan/3 can change the Higgs boson mass by nearly 20 
GeV. This is the most important way in which the parameter tan/3 affects the low-energy 
phenomenology of the model, and it is due to the cos^ 2(3 term in the tree-level part of 
the boundary condition on the Higgs quartic coupling A, eq. ([3]). When Ms gets larger, 
however, the low-energy value of A is dominated by the radiative correction induced by the 
RG evolution, and the effect of the variation in tan/3 becomes less important. 
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Mg (GeV) 



Figure 1: The Higgs boson mass as a function of the scalar mass Ms for = 0, ji = mi/2 = 
500 GeV and two values of tan The dashed lines correspond to a ±lcr variation of the top 
quark mass. 




Figure 2: The Higgs boson mass as a function of the ratio Xt/Ms for Ms = 10 TeV, 
H = mi/2 = 500 GeV and two values of tan/5. The dashed lines correspond to a ±lcr 
variation of the top quark mass. 
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Figure 3: The Higgs boson mass as a function of the scale in various approximations for 
Ms = 10^ and 10^° GeV. For the meaning of the different curves see the text. 

Fig. [2] shows the dependence of on the stop mixing term Xf, for Ms = 10 TeV and 
the other SUSY parameters chosen as in fig. [H It can be seen that, for not too large values 
of Xt/Ms, the dependence of Mh on Xt is milder than in the usual MSSM plots with stop 
masses of the order of 1 TeV (see, e.g., fig. 1.4 of the second paper in ref. [26]). This is due 
to the fact that the threshold correction to the Higgs quartic coupling in eq. ([6]) is computed 
in terms of the running top Yukawa coupling ht{Ms), and the RG evolution up to the large 
scale Ms has the effect of suppressing ht- 

It is also useful to investigate how the one-loop result for Mh depends on the renor- 
malization scale at which we stop the RG evolution of the Lagrangian parameters and 
compute the physical mass of the Higgs boson. In an ideal, all-loop calculation the result for 
the physical mass would not depend at all on the choice of such scale. However, since we are 
truncating our calculation at the one-loop order there will be a residual scale dependence, 
that we can take as a lower bound on (but not necessarily as a full estimate of) the uncer- 
tainty associated with higher-order corrections. In fig. [3] we show the dependence of Mh on 
the scale for the two choices of the heavy scalar mass Ms = 10^ GeV and Ms = 10^° 
GeV (the other relevant parameters are chosen as tan/3 = 30, At = 0, fi{Mz) = mi/2 = 500 
GeV). In each set of curves the dashed line represents the tree-level result, as in eq. ([7]) 
with 6^^ and 6^ set to zero; the dotted line represents the one-loop result in which is 
expressed in terms of the physical top mass Mf, the solid line represents the one- loop result 
in which is expressed in terms of the running top mass mt{Qw) (as is done in figs. [Hand 
[2]). It can be seen that the tree-level result shows a marked dependence on Qw; due to the 
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scale dependence of A in eq. ([7]). A non- negligible dependence is also present in the one-loop 
result based on the physical Mt. On the other hand, the residual scale dependence is very 
small when the running top mass is used in the one-loop correction. This is reminiscent of 
the situation in the MSSM, where the use of the physical top mass in the one-loop corrections 
leads to an excessively high estimate of the light Higgs boson mass, compensated for by large 
negative two-loop corrections, whereas the two-loop corrections are much smaller when the 
running top mass is used in the one-loop part. We also note that the three determinations 
of Mh in fig. [3] are in good agreement with each other for the particular choice Qw ~ ^t- 
The tree-level chargino and neutralino mass matrices read 

/ 



M 



c 



M2 
9dV 



9uV 



M 



N 



Ml 






M2 



9d-" 

9dV 
V2 



_'JdZ 9dJi 



V 



V2 



9dV 

gu V 

V2 







(9) 



The values of the tree-level chargino {X12) neutralino (x? 234) masses depend on the 
renormalization scale at which the various parameters in eq. ([9]) are computed. To reduce this 
scale dependence we include the one-loop corrections to the chargino and neutralino masses, 
once again adapting the MSSM formulae of ref. [14] to the Split-SUSY case (see the appendix 
for details). In figs. H] and [5] we show the masses of charginos and neutralinos, respectively, 
as function of the scale Qw, for Ms = 10^ GeV, tan/3 = 30 and fi{Mz) = rni/2 = 500 GeV. 
In each plot the dashed lines represent the tree-level result and the solid lines represent the 
one- loop result. It can be seen that the inclusion of the radiative correction improves the 
scale dependence of the masses of charginos and neutralinos. For the choice of parameters 
considered in figs. H] and [5] the only exception is the mass of the lightest neutralino, which is 
mostly bino: the tree-level value of m^o has already very little dependence on the scale, and 
the one-loop corrections are negligible. 

Finally, the physical gluino mass Mg is related to the MS parameter M3 by 



Mg 



MsiQ) 



^fl2 + 91ogS 



4tt 



Mi 



(10) 



In the analysis of the constraints on the Split-SUSY parameter space that we present 
in sections [3] and H] we will be mostly interested in phenomena that involve charginos and 
neutralinos. For this reason we find it convenient to choose by default a value of Qw that 
is representative of the masses in the chargino/neutralino sector, i.e. = ^/J^M^. While 
a priori this might not be the best choice of renormalization scale for phenomena involving 
the Higgs boson, fig. [3] shows that the Higgs boson mass is not significantly affected by the 
choice of Qw when we compute the one-loop corrections in terms of the running top mass. 
On the other hand, in many models of SUSY breaking the gluino can be quite heavier than 
the other gauginos, thus we will choose Qw = ^3 when discussing the gluino mass and 
decays. 
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Figure 4: The chargino masses as a function of the scale Qw for Mg = 10^ GeV, tan f3 = 
30, fj, = mi/2 = 500 GeV. The dashed hues are tree-level results and the solid lines are 
one-loop results. 
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Figure 5: Same as fig. |4]for the four neutralino masses. 
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Figure 6: The evolution of the gaugino mass parameters from the GUT scale to the weak 
scale in the universal scenario for Ms = 10^ and 10^° GeV. 

2.2 Patterns of soft SUSY-breaking gaugino masses 

The soft SUSY-breaking gaugino mass parameters Mi^2,3 entering the chargino, neutralino 
and gluino masses are determined via one loop RGEs once their values at the GUT scale are 
fixed. If one assumes universality of these parameters, Mi = M2 = M3 = Tfi\i2 at AfQurj,, as 
is done e.g. in the minimal Supergravity (mSUGRA) model [18], the ratios between the weak 
scale values are simply related to the squares of the gauge couphng constants = gf /{Att) 
by Ml : M2 : M3 = ai : 02 : These values also depend on the intermediate scale Ms below 
which the contributions of the sfermions and of the heavy Higgs doublet are decoupled 
from the RGEs. The evolution of Mi 2,3 is displayed in fig. [6l where a common GUT value 
1^1/2 = 500 GeV has been assumed and two values for the scalar masses, Ms = 10^ GeV 
and Ms = 10^° GeV, have been used. At the scale Mz one obtains the ratios Mi : M2 : M3 = 
1.0 : 2.0 : 7.8 for Ms = 10^ GeV, while one has Mi : M2 : M3 = 1.0 : 2.3 : 10.2 for Ms = 10^° 
GeV. If the scale Ms had been set to 1 TeV, as in the usual MSSM with scalar masses of 
the same order as the gaugino masses, one would have obtained Mi : M2 : M3 ~ 1:2:7. 

Since for heavy scalars the number of basic input parameters of the model is rather 
small, one can relax the assumption of a universal gaugino mass at the GUT scale and still 
have a rather predictive model which, in many instances, could lead to a rather different 
phenomenology. In this paper, rather than performing a general scan with 4 or 5 input 
parameters (i.e. Mi, M2, M3, fi and possibly Ms) we will consider two specific SUSY models 
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in which the boundary conditions for the gaugino masses at the GUT scale are different from 
those of the universal scenario. This will simplify our numerical analysis and, at the same 
time, allow us to address the new interesting phenomenology induced by the non- universality 
of the gaugino masses. 

The first scenario that we consider is a gravity-mediated SUSY-breaking scenario in 
which, to lowest order, the gaugino masses arise from a dimension-5 operator 

£ oc — A A , (11) 

where A"^'^ are the gaugino fields and -F$ the auxiliary component of a left-handed chiral 
superfield $ which couples to the SUSY field strength. In the usual mSUGRA model with 
SU(5) grand unification the SUSY-breaking field F$ is a singlet under the unifying gauge 
group, leading to universal gaugino masses. However, the chiral superfield $ can sit in 
any representation of the symmetric product of the adjoint group [20]. In the case of SU(5) 
symmetry, Fq, could belong to an irreducible representation which results from the symmetric 
product of two adjoint s 

(24 ® 24),y„^etric = 1 © 24 © 75 © 200 (12) 

Once the neutral component of -F$ has acquired a vacuum expectation value, {F<^)ab = 
VaSabi the vevs Va determine the relative magnitude of the soft SUSY-breaking gaugino mass 
parameters Ma at the GUT scale [19]. These are shown in the left-hand side of tableland, 
as can be seen, only in the singlet case 1 are these parameters universal. 

Another set of scenarios that we will consider are four-dimensional string models in which 
the source of SUSY breaking resides predominantly in the moduli sector. In particular, in 
the orbifold Oil scenario in which all chiral fields have modular weight equal to unity, the 
boundary conditions for the gaugino mass parameters at the GUT scale are [21] 

Ma ^ ±10-3(6^ + 5^^)^3/2, a = 1,2,3 (13) 

where m^j^ is the gravitino mass, ha the coefficient of the MSSM beta function for the gauge 
coupling constant Qa-, and 5gs fhe Green-Schwartz mixing parameter which is a negative inte- 
ger number in this case. In these scenarios the scalars, with masses rn^ ^ 10"'^ {—^gs) ^3/2' 
are expected to be much heavier than the gauginos and the pattern of gaugino masses, com- 
pared to the universal case, is driven by the parameter 5gs- For the choice 5gs = ~4 one 
obtains at Mqut the mass pattern Mi : M2 : M3 ~ -y : 5 : 1; see table [H 

Other mass patterns can be obtained by choosing different values of the integer parameter 
5gSi and some of them are in fact similar to those of the non-singlet representation models 
shown in table 1. For instance, for 5gs = ~2, one obtains the ratios Mi : M2 : M3 ~ y : 3 : — 1 
which are close to those of the scenario 200. In the orbifold Ol-type scenario [21], which 
differs from the previous one in the fact that the modular weights are not all equal to unity. 
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Q = Mgut 


Q = Mz [Ms = 10^ GeV] 


Q = Mz [Ms = 10^° GeV] 


1 


1:1:1 


1.0 : 2.0 : 7.8 


1.0 : 2.3 : 10.2 


24 


1 : 3 : -2 


1.0 : 6.3 : -15.2 


1.0 : 6.9 : -19.7 


75 


5 : -3 : -1 


1.0 : -1.2 : -1.5 


1.0 : -1.4 : -2.0 


200 


10 : 2 : 1 


2.4 : 1.0 : 1.9 


2.2 : 1.0 : 2.2 


on 


53/5 : 5 : 1 


1.4 : 1.3 : 1.0 


1.0 : 1.1 : 1.0 



Table 1: The ratios of gaugino mass parameters, Mi: M2: M-i, at the renormahzation scales 
Mgut, Mz (with Ms = 10^ GeV) and again Mz (with Ms = 10^° GeV), for the different 
patterns of soft SUSY breaking. 

one could obtain different gaugino mass ratios. However, in many cases, the pattern is similar 
to that of the scenarios discussed above and, for instance, ratios that are close to those of 
the scenario 75 can be obtained for 6gs = —5. 

The relations between the gaugino masses Mi 2,3 at the scale Mqut and at the weak scale 
Mz are summarized in table [1] for the different scenarios. We have used the one-loop RGEs 
for the couplings and assumed a common scalar mass M5 of either 10'^ GeV or 10^° GeV. 
As can be seen, the pattern of the gaugino mass parameters at Mz (and hence the gluino, 
neutralino and chargino masses) can be quite different from the universal case (scenario 1) 
in which for Ms = 10^ GeV one has approximately Mi : M2 : M3 = cti : 02 : «s — 1:2:8. 

In particular, in the scenario 200 where M2 < Mi, the LSP is wino-like for large values 
of the parameter /x, implying that Xi xf ^ire nearly degenerate in mass. Again for large 
H values, the neutralinos Xi a-^id X2 ^^e charginos xf masses are very close to each other 
in the scenario 75 since |Mi| ~ IM2I, while in the scenario 24 the mass splitting between the 
LSP and the states X2 xf can be very large since M2 ~ 6M1. Finally, in the Oil model 
one has M3 < Mi, M2 and the gluino tends to be the LSP unless fi is very small, in which 
case Xi; X2 ^^"^ Xi are higgsino-like and almost degenerate in mass. Note that, in general, 
the weak-scale ratios among the three gaugino masses show a dependence on the scale Ms 
at which the sfermions are integrated out. 

There are many other SUSY models that lead to non-universal gaugino masses, and a 
review has recently been given in ref. [27]. However, in most cases one is very close in 
practice to the patterns that have been introduced above, and the phenomenology of the 
gaugino sector is quite similar. For instance, the weak-scale pattern of gaugino masses that 
emerges from anomaly mediated SUSY-breaking (AMSB) models [28], Mi : M2 ~ 3:1, is 
similar to that of the model 200. In mirage gaugino mediation [29], where SUSY breaking 
is realized in higher-dimensional brane models (a scheme that is realized naturally in the 



13 



so-called KKLT-type moduli stabilization models [30]), one has [27] Mi:M2~l:1.3 which, 
again, is similar to the pattern of the scenario 75. Thus, we believe that the scenarios 
discussed above, with the patterns of gaugino masses of table 1, are representative of a wide 
spectrum of non-universal models. 

2.3 Implementation of the MSSM with heavy scalars in SuSpect 

We have implemented this MSSM scenario with heavy scalars into the RGE code SuSpect 
[11]. This model can be chosen by selecting at the very beginning of the input file suspect2 . in 
or alternatively suspect2_lha. in, the option SHeavy by putting ichoice(l) = 3. The only 
two sets of basic input parameters needed to be set are: 

- The SM basic input parameters, i.e. the electromagnetic, strong and weak couplings, 
the Z boson mass and the third-family fermion masses: 

a{Mz), as{Mz), Gp, Mz, Mt, mb{mb), Mr 

- The additional input parameters specific for this model: 

fi{Mz), Mi(Mgut), M2(Mgut), Ms{Mgvt), Ms, tanp{Ms), At{Ms). 

All the other parameters of these files are irrelevant. 

The routine SHeavy . f performs the RG evolution for the gauge couplings, the third gen- 
eration fermion Yukawa couplings, the gaugino-higgsino-Higgs boson couplings {gu,d, gu,d), 
the gaugino mass parameters, the /i parameter, and the quartic Higgs coupling A. This 
routine also contains all the relevant one-loop radiative corrections. 

The output file contains the physical chargino and neutralino masses and the elements 
of the mixing matrices U, V and (see the appendix for details), as well as the physical 
masses of the gluino and of the lightest Higgs boson. They are computed at a scale Qw 
set by default to Qw = y^JiM^, except for the gluino mass which is computed at the scale 
M3. The masses of all the sfermions and of the heavy Higgs bosons are considered to be 
degenerate and set to M5, while the various mixing angles (in the third-generation-sfermion 
and Higgs sectors) are set to zero. 

Note that, for the phenomenological analyses that will be presented in the next two 
sections, we also needed to adapt to the Split-SUSY scenario the two programs HDECAY 
[16] and SDECAY [17], which compute the decay widths and branching ratios of the MSSM 
Higgs bosons and of the SUSY particles, respectively. These programs use the output given 
by SuSpect for the soft SUSY-breaking parameters, the mixing matrix elements and the 
sparticle and Higgs masses, but they calculate internally the various couplings. In particular, 
the Higgs couplings to neutralinos and charginos (which are different from the usual MSSM 
case) are hard-coded and need to be adapted. 
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3 Collider and Dark Matter constraints 



In this section we analyze the constraints on the MSSM with heavy scalars, first from collider 
searches and high-precision data [22] and then from cosmological data, in particular the relic 
density measurement of DM by the WMAP satellite [31] and the gluino lifetime [32-34]. 
A special attention will be given to the non-universal gaugino mass scenarios discussed in 
the previous section, where several new features compared to the universal case appear. 
The DM and some collider constraints for non-universal gaugino masses have been discussed 
in ref. [35] in the usual MSSM with light scalars, while some implications of DM in the 
Split-SUSY scenario have been studied in Refs. [15,36]. 

3.1 Constraints from collider data 

In the scenario in which a universal gaugino mass mi/2 is assumed at the GUT scale, leading 
to the approximate relation Mi : M2 : M3 ~ 1 : 2 : 8 at the weak scale, the strongest 
experimental bound is due to the negative search of charginos at LEP2 up to energies of 
^/s ~ 208 GeV. From pair production of the lightest chargino, e^e~ XiXii obtains 
the mass bound [22] 

m-± > 103 GeV . (14) 

The chargino mass bound in eq. f|T^ is valid only if the mass splitting between the lightest 
chargino and neutralino is large enough, AM = m+ — m.yO > a few GeV. For small AM 
values, as is the case when the LSP neutralino is higgsino- or wino-like, m^± ^ m^o ^ or 
M2, the bound becomes weaker and, for very heavy scalar fermions, one has m-± > 92 GeV 
[22]. Using the weak-scale gaugino mass relations above, the bound on the lightest chargino 
mass translates into a lower bound on the LSP mass, m^o > 50 GeV if the LSP neutralino 
is bino-like, in which case one has m^o ~ Mi ~ — ^n^^y^f] Xi is thus too heavy to be 
kinematically accessible in invisible Z boson decays at LEPl, Z —>■ XiXi- the higgsino- and 
wino-like regions for the LSP neutralino the mass bound is higher, m^o > 92 GeV, as Xi is 
almost degenerate with the chargino xf discussed above. Furthermore, the constraint on 
M2 from the bound on m+ translates into a constraint on the gluino mass, nig ~ M3 > 350 
GeV, which is higher than the direct bound from Tevatron searches when scalar quarks are 
very heavy [22] 

m-g > 195 GeV . (15) 

The constraints on the [M2,/i] parameter space from LEP2 and Tevatron negative searches 
of SUSY particles are summarized in fig. [71 Here and in the following M2 and n have to 
be interpreted as MS parameters evaluated at the renormalization scale Q = Qyf. Besides 
the process Z — > XiXi^ which contributes to the invisible Z boson decay that is tightly 
constrained by LEPl data, P^^^ < 2 MeV [22], we have imposed that the cross sections for 
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Figure 7: Constraints on the [M2,/i] parameter space for the scenario 1 with universal 
gaugino masses at the GUT scale and with Ms = 10^ GeV, At = and tan/3 = 30; the 
different colors correspond to the regions excluded by the LEP bounds on the partial decay 
width r(Z — > XiXi) on the production cross sections e~^e~ — > xfxi^ X1X2? XiXs- 

the three processes e~^e~ xfXi^ e~^e~ — X1X2 — X1X3 ^^e smaller than 50 fb 

which, given the collected luminosity of £ ~ 100 pb~^ at the cm. energy ^/s ~ 208 GeV, 
corresponds to less than 5 events. As can be seen, the strongest constraint is by far due to 
eq. ( fT^ . i.e. the green (light grey) area. The process e"'"e~ — > X1X2 where the LSP and the 
next-to-lightest neutralino are pair produced, although more favored by phase space, does 
not add much information as the cross section is generally much smaller. In the small blue 
(dark grey) oval region on the left of the plot X2 is mostly wino, while in the region with 
large M2 and small /i on the right of the plot X2 is mostly higgsino. 

All these bounds can be adapted to scenarios in which the boundary conditions for the 
gaugino masses at the high scale are different. However, in some cases, their impact can be 
widely different compared to the universal scenario, as shown in fig. [8] where the constrained 
[M2,/i] parameter space with Ms = 10'^ GeV is displayed for the four scenarios 24, 75, 200 
and Oil. We have assumed /i > but a similar pattern is obtained for /i < 0. 

In the non-universal scenario 24 one has M2 : Mi ~ 6 : 1 for the wino and bino masses 
at the weak scale and, as in the scenario 1, the chargino mass bound in eq. (fT^ leads to 
the strongest constraint and rules out the entire fi, M2 > 100 GeV range. For large values 
of the wino mass, M2 > 300 GeV, and small values of fi, the phase space for the process 
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e+e~ — s> X1X3 is still open and a small additional region is ruled out. The process where 
the second neutralino is produced in association with the LSP, e^e~ — > X1X2 5 ^^so plays a 
role at even larger values of M2 and rules out another portion of the parameter space. If 
the lighter chargino and neutralinos are gaugino-like, M2, Mi <^ the LSP neutralino is 
bino-like and the chargino mass bound in eq. (fT4|) translates into the relatively weak bound 
m^o > 17 GeV. This opens the possibility that the decay mode Z — >• XiXi contributes to the 
invisible decay of the Z boson. For moderate values of /i, for which the LSP has a higgsino 
component and hence sizable couplings to the Z boson, the constraint Pinv ^ 2 MeV rules 
out a region in the [M2, ^] plane that is not excluded by any other process. 

In the scenario 75 the masses of the lightest chargino and of the two lightest neutralinos 
are very close to each other, since Mi : IM2I ~ 1 : 1.2. Therefore, chargino pair production, 
which has the largest cross section among the three processes with a similar phase space 
e^e~ — > xfxf^ X1X2 X2X21 leads to the strongest constraint on the [M2, n] parameter 
space. However, since the weak-scale value of M3 is also close to the values of Mi and M2, 
the constraint rrig > 195 GeV from negative searches of gluinos at the Tevatron plays an 
important role, in contrast to the scenarios 1 and 24 in which M3 > 3-4 M2 and the bound 
of eq. ( ITSj) is superseded by that of eq. ( fT4l) . Since M3 ~ 1.3 M2, the entire area M2 > 150 
GeV is ruled out by eq. (fT5|) independently of the value of fi. 

A similar pattern occurs in the scenario 200 where the mass parameters Mi,M2 and 
M3 are also close to each other. However, since here the ratio M3/M2 is larger than in the 
scenario 75, the chargino mass bound is stronger than the bound from gluino searches. Note 
also that, in this scenario, the LSP is wino-like for large values of n and the chargino mass 
bound translates to m^o > 100 GeV for the LSP and the two particles are almost degenerate 
in mass (this also holds in the higgsino-like region for the LSP). 

Finally, the string-inspired scenario Oil is interesting in many respects. As in the uni- 
versal scenario 1 chargino pair production is the most constraining of all LEP production 
processes and rules out the same area of the [M2,/i] parameter space. However, since here 
the smallest of the gaugino masses is the gluino mass M3, the Tevatron bound in eq. (|T5|) 
rules out a significant portion of the parameter space, namely M2 < 230 GeV. Furthermore, 
the gluino tends to be the LSP for large values of fi when the lightest neutralino is gaugino- 
like. A strongly interacting stable particle is disfavored by cosmological data, therefore this 
area of the parameter space has to be excluded. 

Another important collider constraint comes from the negative search of Higgs bosons at 
LEP2 where a lower bound, Mh ^114 GeV, has been set on the mass of a SM-like Higgs 
boson [22]. As previously discussed, the Higgs boson in the MSSM with heavy scalars is 
SM-like but its mass is generally larger than 130 GeV (for M5 > 10^ GeV and tan/3 large 
enough), therefore this constraint does not apply in our case. Note, however, that there is a 
chance that such a Higgs particle could be observed at the Tevatron, either in the processes 
qq HW bUu for Mh ~ 130 GeV or in the process gg ^ H ^ WW^*'^ for Mh ~ 160 
GeV, if a large integrated luminosity is collected; see for instance ref. [26] for details. 
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Figure 8: The same as in fig. [7] but for tlie non-universal scenarios. Where relevant, 
effect of the gluino mass bound from Tevatron searches is also shown. 
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Finally, we summarize the constraints from high-precision data and rare decays. Be- 
cause the sfermions are very heavy, the SUSY-particle contributions to the anomalous mag- 
netic moment [g — 2)^ of the muon (which occur essentially through smuon-neutralino and 
sneutrino-chargino loops, with possible very small contributions from Higgs bosons) and to 
the rare decay of the 6-quark into a strange quark and a photon BR(6 57) (which at 
leading order occur via loops involving the charged Higgs bosons and top quarks as well as 
charginos and top squarks) are extremely small, and one is left only with the SM contribu- 
tions. The effect of SUSY particles on the high-precision electroweak observables measured 
at LEP, SLC, Tevatron and elsewhere is also very tiny in the case of heavy scalars. Indeed, 
the dominant contributions to these observables, in particular the W boson mass and the 
effective weak mixing angle siii^ 9\y, enter via a deviation from unity of the p parameter 
(which measures the relative strength of the neutral to charged current processes at zero 
momentum transfer). The sfermions and the non-SM Higgs bosons are heavy enough that 
their contribution to this quantity is negligible. The chargino and neutralino contributions 
are small because the only terms in the mass matrices that could break the custodial SU(2) 
symmetry are proportional to Mw- This has been verified explicitly in the case of universal 
gaugino masses in ref. [37], where it has been shown that only when charginos and neu- 
tralinos have masses very close to the experimental lower bounds can they affect, and only 
slightly, the electroweak observables. 

3.2 The Dark Matter constraint 

As deduced from the WMAP satellite measurement of the temperature anisotropies in the 
Cosmic Microwave Background, cold Dark Matter makes up approximately 25% of the energy 
of the Universe [31]. The DM cosmological density is precisely measured to be 

Qj^uh' = 0.11i;°;°08 , (16) 

which leads to 0.088 < QDuh"^ < 0.128 at the 3a level. The accuracy is expected to be 
improved to the percent level by future measurements at Planck. 

As is well known, the LSP neutralino is an ideal candidate for the weakly interacting 
massive particle that is expected to form this cold DM [5] and in some areas of the SUSY 
parameter space the Xi cosmological relic density, which is inversely proportional to the 
neutralino annihilation cross section (Jann = cr(XiXi ~^ SM particles), falls in the range 
required by WMAP. In the MSSM with heavy scalars, there are essentially only three regions 
(see later for a fourth possibility) in which this constraint is satisfied: 

- The "mixed region" in which the LSP is a higgsino-gaugino mixture. Mi ~ which 
enhances (but not too much) its annihilation cross sections into final states containing 
gauge and/or Higgs bosons and top quarks, XiXi ~* W~^W~ , ZZ, HZ, HH and ti. 
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- The "pure higgsino" and "pure wino" regions, in which the LSP is almost (but not ex- 
actly) degenerate in mass with the lightest chargino and the next-to-lightest neutralino, 
leading to and enhanced destruction of sparticles since the xt ^ X2 co-annihilation cross 
sections are much larger than that of the LSP; this solution generally requires LSP 
masses beyond 1 TeV. 

- The "if-pole" region in which the LSP is rather light, m^o ~ ^Mh, and the s-channel 
H exchange is nearly resonant allowing the neutralinos to annihilate efficiently. 

The latter scenario has been discussed in the usual MSSM [38] in which the Higgs boson, 
which has a mass below ~ 130 GeV, decays mostly into bb pairs. However, if the common 
scalar mass Ms is very large, as is the case in the present discussion, the Higgs mass will 
be larger reaching, for instance, the value Mh ~ 150 GeV for Ms ~ 10^*^ GeV. In this case 
the dominant decay mode of the Higgs boson is a pair of W bosons, one being on mass-shell 
while the other is virtual and decays into two massless fermions, H — s> WW* — > Wff. 

The partial decay width for this decay is given by [39] 



with X = M^/Mfj. Using the two-body-like formula given above and the approximation of 
ignoring the kinematical effects induced by the three-body final-state phase space, we have 
implemented this additional channel in a DM Fortran code based on ref. [40] which is linked 
to the program SuSpect. Using this routine we have scanned the [Mo, yu] parameter space 
and determined the regions in which the WMAP constraint is fulfilleclj. This is performed 
not only in the scenario with a universal gaugino mass mi/2 at the GUT scale but also in 
scenarios with different boundary conditions. 

In fig. [9] we display the area in the [M2, fj] parameter space in which the WMAP constraint 
is fulfilled; a common scalar mass Ms = 10^ GeV is chosen and a universal gaugino mass at 
Mgut is assumed; here and in the subsequent discussions we will fix the value of tan/3 to 
tan/3 = 30. The green (light grey) area in the left and bottom parts of the figure denotes 
the region excluded by the collider data discussed in the previous subsection. The peak for 
small M2 values, M2 ~ 2Mi ~ Mh, is due to the s-channel exchange of the Higgs boson, 
XiXi ~^ H . For the mass value obtained here, Mh ~ 130 GeV, the Higgs boson mainly 
decays into bb final states while the H —>■ WW* channel, which has also been included, has a 
smaller branching ratio and does not play a leading role. Between the two bands of the peak 

^For some selected points of the parameter space, we have verified that the results that we obtain for 
ilh^ using the routine of ref. [40] are in a relatively good agreement with those obtained with the program 
micrOMEGAs [41]. 
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Figure 9: The regions of the [M2,fi] parameter space in which the WMAP constraint is 
fulfilled (red area) for a common scalar mass value Ms = 10^ GeV, tan/? = 30 and the 
assumption of a universal gaugino mass at the GUT scale. The green area on the left and 
the bottom is the one excluded by direct searches of SUSY particles. 

one is too close to the Higgs mass pole, and the LSP annihilation is too efficient leading to 
a too small Qh'^. The peak reaches up to /i ~ 600 GeV, a value beyond which the LSP is 
almost bino-like and its coupling to the Higgs boson is too small (the Higgs prefers to couple 
to a higgsino-gaugino mixture) to generate a sizable annihilation cross section. 

For larger fi and M2 values there is an almost straight band in which fi ~ Mi and the 
LSP is a bino-higgsino mixture with sizable couplings to W, Z and Higgs bosons, allowing 
for reasonably large rates for neutralino annihilation into XiXi ~^ W~^W~, ZZ, HZ and HH 
final states. For instance, for M2 ~ 300 GeV and fi = 200 GeV, the annihilation cross 
section is mostly due to the WW and HH final states (~ 40% for both channels) and, to a 
lesser extent, the ZZ and ZH final states (~ 10% each). For slightly larger /x and M2 values 
there is a jump due to the opening of the XiXi ~^ tt channel, which then dominates the 
annihilation cross section. Above the band and below the band, the LSP couplings to the 
various final states are either too strong or too weak to generate the relevant relic density. 
For fi values close to 1 TeV and even larger values of M2 there is a wider area in which the 
WMAP constraint is also fulfilled. In this region the LSP is almost a pure higgsino and a 
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Figure 10: The same as in fig. [9] but with Ms = 10^° GeV (left) or /i < (right) 

correct fi/i^ can also be obtained thanks to the co-annihilation of the LSP with the xf ^iid 
X2 states. For lower /i values and M2 still very large the LSP co-annihilation with xf and 
X2 is too strong and leads to a too small Qh"^. 

Figure [10] is similar to the previous one, with the difference that Ms = 10^° GeV (left 
pane) or the sign of the parameter /i is reversed (right pane). The figure in the left-hand 
side shows similar features as fig. [HI except that the peak due to the s-channel Higgs boson 
exchange is shifted to a slightly higher M2 value, M2 ~ 2Mi ~ M^ ~ 150 GeV. In this 
case the new annihilation channel XiXi H ^ WW* — > Wff discussed before gives a 
significant contribution to the total cross-section. For negative fi values the [M2,/i] area of 
the parameter space that leads to the desired Qh"^ is also similar to fig. [H The only difference 
is that the Higgs peak reaches up to ~ 400 GeV only, as for negative fi values the LSP 
becomes bino-like more quickly than in the positive case and its couplings to the Higgs boson 
are thus smaller. 

Figure [TT] shows the [M2,/i] area that is compatible with WMAP results for Ms = 10^ 
GeV in the case where the boundary conditions for the gaugino masses at the GUT scale 
are not universal. In the scenario 24 the same trend as for the universal case occurs but 
with two major differences. First, the Higgs peak is now shifted to M2 ~ 400 GeV, a 
mere consequence of the fact that the ratio of the weak-scale wino-to-bino masses is much 
larger in this model, M2 : Mi ~ 6, than in the universal case, M2 : Mi ~ 2; for the same 
reason, the band at large M2 values is shifted downward compared to fig. [9] The other 
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major difference with the universal scenario is that, despite the constraint on the invisible Z 
decay width YiZ XiXi) ~ 2 MeV and the one from chargino and neutralino production 
at LEP2, the possibility that the LSP mass is close to \Mz is not excluded. In this case, 
the LSP annihilation channel XiXi ~^ ^ ^ ff can become resonant, thus generating the 
required cosmological relic density. This leads to a peak similar to the one due to Higgs 
boson exchange but at smaller M2 values, M2 — 6 Mi ~ 3 Mz ~ 300 GeV. 

In the scenario 75 one has Mi : M2 : M3 ~ 1 : —1.2 : —1.5 at the weak scale, so that the LSP 
is in general close in mass to the lightest chargino and the next-to-lightest neutralino, and 
co-annihilation of these states plays a very important role. Indeed, in the thin straight line 
below /i, M2 < 1 TeV the LSP is dominantly bino-like (e.g. ~ 85% for = M2 ~ 600 GeV) 
with a small higgsino component, and the WMAP Qh^ range is obtained with the efficient 
co-annihilation of X^Xi'jXi'Xi' (~ 50%) and XiXi" 5X1X2 (~ 25%) with the remaining part 
due to XiXi annihilation. In the large band with fi = 900-1100 GeV for M2 ^1.2 TeV 
the LSP is very close to be a pure higgsino state (e.g. ~ 99% for M2 ~ 2/^ ~ 2 TeV), but 
the co-annihilation cross sections are still reasonably small, the mass difference between the 
LSP and the other chargino/neutralino states being large enough. All other areas, including 
the Higgs peak and the mixed gaugino-higgsino areas that appear in the universal case, are 
ruled out because of a too efficient co-annihilation of the LSP. 

The same situation occurs in the scenario 200, in which the weak-scale gaugino mass 
parameters are given by Mi : M2 : M3 ~ 2.4 : 1 : 1.9 leading to Xi ^^^d xf states which are 
wino- or higgsino-like and almost mass-degenerate in most of the parameter space. In this 
case, co-annihilation with xf and X2 too efficient except for the band with fx = 900-1100 
GeV and M2 > 1.2 TeV as in the scenario 75. In this case the band is extended at larger 
H values by a strip in which the LSP is a wino-higgsino mixture, but, again, the bulk of 
the relic cosmological density Qh"^ is generated through co-annihilation of the LSP with the 
lightest chargino Xi"- 

Finally, in the case of the Oil string model, which leads to a weak-scale gluino mass 
parameter that is smaller than the wino and bino mass parameters. Mi : M2 : M3 ~ 1.4 : 1.3 : 1 
for Ms = 10^ GeV, the gluino is the LSP in a large part of the parameter space and this 
blue (dark grey) area should therefore be excluded. For M2 smaller than 1 TeV the LSP is 
higgsino-like at low /i values and co-annihilation of the LSP with xf and X2 ^^o efficient, 
while at high /i values the LSP is the gluino. Only in a relatively narrow band, /i ~ 800- 
1100 GeV, similar to the ones observed in the scenarios 75 and 200, does one obtain the 
Qh"^ range required by WMAP, with the LSP being an almost pure higgsino state with a 
relatively large mass; for instance, for M2 ~ 2.5 fi ~ 2.5 TeV, one has m^o ~ m^o ~ m^± ~ 1 
TeV and Qh"^ is almost exclusively generated by co-annihilation (95%). A very interesting 
feature occurs at the border between the gluino-LSP area and the band with fj, ~ 0.8-1 TeV: 
the main channel that leads to the required range for the relic density is gluino annihilation 
into gluon and quark pairs, gg — > qq, gg (the co-annihilation gx qq is suppressed as it 
is mediated by the super- heavy squarks). For instance, in the point M2 ^1.5 TeV and 
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Figure 11: The regions of the [M2,yU.] parameter space in which the WMAP constraint is 
fulfilled for Ms = 10^ GeV in the various models with non-universal GUT-scale boundary 
conditions for the gaugino mass parameters. 

/i ~ 900 GeV one obtains the value flh"^ ~ 0.1 with the gg gg (45%) and gg — > qq (35%) 
reactions. This is one of the rare examples within constrained SUSY models where gluino 
co-annihilation is at work. 
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3.3 The gluino lifetime 

In this section we summarize the constraints that can be obtained on the parameter space, 
and in particular on M5, from the requirement that the gluino lifetime does not exceed 
the age of the Universe, Tg < 14 Gyear, a possibility that is excluded from the absence 
of anomalous isotopes. We extend the discussions of Refs. [32-34], held in the context of 
universal gaugino masses, to the various scenarios with non-universal boundary conditions 
at Mgut- 

The total decay width of the gluino, Tg = h/Tg, has been calculated following ref. [34], 
where the large logarithmic corrections that appear for very heavy scalars and are controlled 
by the strong coupling and by the top Yukawa coupling are resummed with an effective 
Lagrangian approach. Both the three-body decays into charginos or neutralinos and a quark- 
antiquark pair through the exchange of heavy squarks, g — > qq* qqXii ^-nd the loop induced 
decay into a neutralino and a gluon, g — > gx^i have been included (see section 14.31 for more 
details). The gluino lifetime is approximately given by [34] 



where is a normalization factor which is generally of order unity if phase-space effects 
are ignored. This equation exhibits the main trend: the gluino lifetime is larger for higher 
values of the scalar mass Ms and smaller values of the gluino mass rrtg. In Refs. [32-34] it 
has been shown in the universal scenario that for = 0(10^3) GeV the gluino is almost 
stable as its lifetime is larger than the age of the Universe. 

The gluino lifetime Tg is displayed in fig. [12] as a function of Ms for various values of 
the gluino mass in the two non-universal scenarios 24 and Oil; for each value of and 
hence of M2 (the latter obtained via RG evolution from the boundary conditions specific to 
the scenario) the relevant value of yU. is obtained by requiring that the LSP relic abundance 
Vt~)Ji^ falls in the range allowed by WMAP (see fig. [TTl) . As can be seen, in the scenario 24 
a gluino lifetime of the order of the age of the Universe, Tg = 14 Gyear (the horizontal line), 
is obtained for Ms = 10^^ to 10^^ GeV for rrig = 1-3 TeV. The results are thus similar to 
the universal scenario discussed in Refs. [8,34]. In fact, the same results are also obtained 
in the non universal scenarios 75 and 200 and we refrain from showing them again. 

A slightly different situation occurs in the scenario Oil, in which the correct DM relic 
density can be generated when the gluino is close in mass to the higgsino-like neutralino LSP. 
In this case the only allowed decays of the gluino are the three-body decays into two light 
quarks and the higgsino-like chargino and neutralinos, and the loop-induced two-body decay 
into the higgsino-like neutralinos and a gluon (see section 4.3). These decays are suppressed 
by phase space: F cx: {rrig — m^Y fo^ the three-body decays and F oc {rrig — m^)'^ for the 
radiative decay. Therefore, smaller values of Ms are required for the gluino lifetime to be 
compatible with the age of the Universe. This is exemplified in the right-hand side of fig. [T2|, 
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Figure 12: The gluino lifetime (in seconds) as a function of the common scalar mass Ms for 
the two non-universal scenarios 24 and Oil for various values of nig and the parameter 
fixed from the DM requirement; the horizontal line is for the upper bound of 14 Gyear. 

where one can see that for rrig = 0{1 TeV) (with /i ~ 1 TeV for Q^h"^ to fall in the WMAP 
range) a gluino lifetime of Tg = 14 Gyear is already reached for Ms < 10^^ GeV. 
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4 Decays of the Higgs and SUSY particles 



4.1 Higgs boson decays 

As discussed in section 2.1, in the MSSM with heavy scalars the Higgs boson is SM-hke, 
except that its mass is constrained to be in the range Mh ~ 130-160 GeV for Ms = 10^- 
10^^ GeV and tan/3 = 30 (as will be assumed throughout this section). It will thus decay 
mostly like the SM Higgs particle in this mass range [26,42]. For not too large Ms values 
for which Mh ^ 130 GeV, the Higgs boson decays into a large variety of channels, the main 
modes being by far the decay into bb pairs with a branching ratio of ~ 90% followed by 
the decays into cc and t~^t~ pairs with branching ratios of ~ 5%. Also of significance, 
the top-loop mediated Higgs decay into gluons which for Mh around 130 GeV occurs at 
the level of few percent. The top- and 14^- loop mediated 77 and Z7 decay modes are very 
rare, with branching ratios of O{10"^). However, these decays lead to clear signals and 
are theoretically interesting, being sensitive to new electrically charged particles such as 
charginos. For values of Ms large enough that Mh ^ 140 GeV the Higgs bosons decay into 
WW and to a lesser extent ZZ pairs, with one of the gauge bosons being virtual below the 
threshold. For Ms > 10^° GeV, which leads to Mh ^ 150 GeV, the Higgs boson decays 
almost exclusively into two real W bosons; the decay H ZZ* is strongly suppressed as 
one of the Z boson must be virtual. In all cases, the Higgs boson is very narrow, as its total 
decay width does not reach the 1 GeV level. 

There are, however, two situations in which the Higgs boson might have decays that are 
slightly different from those of the SM Higgs particle. First, for very light LSP neutralinos, 
the invisible decay H — * XiXi [43] might be kinematically accessible. In the scenario with 
universal gaugino masses leading to M2 ~ 2 Mi at the weak scale, this decay occurs for values 
of M2 small enough to have the phase space needed for the decay to occur, Mh ^ 2m^o ~ 
2 Ml. In the left-hand side of fig. [13] the areas in the [M2,/i] parameter space in which the 
branching ratio BR(i7 XiXi) is larger than 1%, 5% and 10% are shown for Ms = 10"^ 
GeV and tan/3 = 30; the area in which the neutralino relic density is in the range required 
by WMAP is also displayed, as well as the area excluded by collider bounds. For M2 values 
in the range 120-150 GeV and small fi values, the branching ratio is of the order of 5% 
and sometimes 10% and is therefore measurable at the ILC. The branching ratio drops with 
increasing n values since for /x ^ M2 the Higgs-LSP coupling becomes too small. In the 
non-universal scenario 24 with the weak-scale relation M2 ~ 6M1 between the wino and bino 
masses, fig. [I3] (right), the corresponding areas in the [M2, /i] parameter space are larger as a 
result of a larger phase space allowed for the invisible Higgs decays. In particular, branching 
ratios larger than 10% are possible in a significant portion of the parameter in which the DM 
constraint is also fulfilled. Thus, if by chance it is the H or Z boson pole which provides the 
correct value of Qh"^, the invisible decay branching ratio BR(if —>■ XiXi) could be measured 
at the ILC and would allow to access directly to the Higgs-LSP couplings. 
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Figure 13: The areas in the [M2,yu] plane in which the invisible Higgs branching ratios 
BR(i? — > XiXi) ill larger than 1%, 5% and 10% in the scenarios 1 and 24 with Ms = 10^ 
GeV and tan f3 = 30; the region in which the relic density is compatible with WMAP is also 
displayed (red region). The program HDECAY [16], adapted to deal with heavy scalars, has 
been used. 



The decays of the Higgs boson into the heavier neutralinos and the charginos are in 
general kinematically closed if one takes into account the LEP2 bounds on the masses of 
these particles. There is however one possible exception: in the scenario 24, because the 
absolute lower limit on the LSP mass is only m^o > 17 GeV, the possibility of the decay 
H X1X2 is s^iii open. This occurs for values of Mi and M2 that lie very close to those 
ruled out by the experimental constraints and, for such values, the requirement that the LSP 
provides the correct relic density is not fulfilled. 

Another possibility for a deviation from the SM predictions for the Higgs decays is when 
the charginos are light enough to contribute to the loop decay i7 — * 77 (contributions to the 
H Z'-f decays are in general much smaller) [44]. Indeed, besides top quark and W boson 
loops, one has to include also the contribution of the chargino loops. However, in contrast 
to the SM particles, the charginos do not couple to the Higgs boson proportionally to their 
masses and the amplitudes are damped by inverse powers of the masses. The chargino 
contributions are thus sizable only for relatively small masses; see also ref. [45] in which this 
topic has been discussed in the Split-SUSY scenario. 

This is exemplified in fig. [HI where the regions in the [M2, /i] plane in which the deviation 
of r{H 77) from SM prediction is larger than 1% and 2% are displayed in the universal 
gaugino mass scenario with M5 = 10^ GeV; the area where the WMAP DM constraint in 
fulfilled has been superimposed. Only for small M2 values, and thus rather light chargino 
states, does the deviation from the SM prediction exceed the level of 2%, which makes it 
potentially observable at the 77 option of the ILC (where one expects the H'-ff coupling to be 
measured at the two-percent level). The sign of the contribution is controlled by the sign of 
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Figure 14: The areas of the [M2, n\ plane in which the deviation of T{H 77) from the SM 
prediction is larger than 1% and 2% for Ms = 10"^ GeV; the area where the WMAP DM 
constraint in fulfilled is also shown. The program HDECAY [16], adapted to deal with heavy 
scalars, has been used. 



/X, thus for yU < one would have a negative shift in r{H 77). The chargino contribution 
drops for higher /i values as, in this case, the chargino xf which gives the most important 
contribution becomes more wino-like and has a weaker coupling to the Higgs boson. 

Since the chargino masses depend only on M2 and fi, the same figure holds for the non- 
universal scenarios; the only difference is that the WMAP-allowed areas for the LSP relic 
density, which have been given in fig. [TTl would be different. In fact, one can see from 
these figures that in the scenarios 75, 200 and Oil the DM constraint is fulfilled only for 
large fj,, M2 values, thus for charginos too heavy to contribute to the 77 decay of the Higgs 
boson. Furthermore, in these scenarios neutralino LSPs that are compatible with WMAP 
data are in general also too heavy to allow for the occurrence of the H XiXi invisible 
decays discussed earlier. 

Note that for very large Ms values, Ms ^ 10^° GeV when the Higgs mass Mh ~ 160 
GeV becomes close to the WW threshold, the branching ratios for both the invisible decay 
and the 77 decay of the Higgs boson become smaller than for Ms = 10^ GeV, as these decays 
have to compete with the H — > WW* channel, which has a much larger decay rate than the 
H ^ bb channel. 
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4.2 Chargino and neutralino decays 

In most cases, the charginos and the neutrahnos (except for the LSP) will decay into lighter 
X states and V = W/Z gauge bosons which can possibly be virtual and subsequently decay 
into two massless fermions, Xi ~^ Xj'^^*^ ~^ Xiff [46]. As the scalar fermions are very heavy, 
their virtual exchange Xi ~^ f f* ~^ ffXj is strongly suppressed and they do not participate 
in the decay processes. The branching ratio for the full final states will thus essentially 
follow that of the gauge bosons, therefore the branching ratios into I = leptons are 
rather small: 20% for the charged and 6% for the neutral decays. There are however two 
interesting features which might occur and which will be briefly discussed below: i) decays of 
the heavier charginos/neutralinos into lighter ones and the Higgs boson, Xi ~^ XjH [47] and 
a) the loop-induced decay of some neutralinos into the LSP and a photon, x° — >• Xil [48]. 
The former process would allow to access directly the Hxx couplings and the latter has an 
interesting experimental signature, a monochromatic photon. 

Adapting the program SDECAY [17] to the case of heavy scalars, we have calculated 
the branching ratios for the decays of the charginos and heavier neutralinos into lighter 
states and the Higgs boson. The result is illustrated in fig. [15], where the branching ratios 
BR(xi — > XjH) are shown as a function of in the two scenarios 1 (top) and 24 (bottom) 
for Ms = 10^ GeV; for a given /i, the value of M2 is chosen in such a way that we sit on 
the left band of the Higgs peak in the WMAP-allowed region of the [M2, /i] plane (see figs. [9] 
and fTTl) . Two main ingredients control the size of the Xi ~^ XjH branching fractions: 

i) The mass difference between the initial and final x states and hence the importance 
of the phase space; in fact, the Xi~Xj mass difference needs to be larger than Mh for the 
Higgs boson to be on-shell, as the H* — > // virtuality would be strongly suppressed by the 
small Hff couplings. 

ii) The initial and final neutralino and chargino states should have different textures as 
to maximize the HxiXj coupling. 

As can be seen from the figures, some decays such as X2 ~^ XiH and X4 ~^ X2H can 
reach the 20% level. In the scenario 1 the other decay modes involving the Higgs boson 
are below the 10% level as they are suppressed either by phase space or by the smaller 
Hxx couplings. In the scenario 24 an interesting feature occurs for intermediate /i values, 
fi = 300-450 GeV, where all charginos and neutralinos are mixed gaugino-higgsino states 
and have masses of the same order. In this range, the dominant decay channels mentioned 
above are kinematically closed, allowing for the phase-space favored decays X34 ^ XiH to 
dominate. 

For the radiative and loop-induced decay Xi ~^ Xj7 (which is of higher order in per- 
turbation theory and thus suppressed by additional powers of the electroweak coupling) to 
occur with a substantial rate, the standard decay modes X^^y x'jH need to be strongly 

suppressed. This occurs when the Xi~X^ mass difference is smaller than Mz and thus Mh, 
so that the decay is a three-body process with a partial width that is suppressed by the 
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Figure 15: The branching ratios for the decays of the heavier charginos and neutrahnos into 
hghter ones and Higgs bosons, BR(xj XjH) as a function of /i in the scenarios 1 (top) 
and 24 (bottom); Ms = 10^ GeV is assumed and M2 is chosen such that the rehc density is 
compatible with WMAP. The program SDECAY [17], adapted to deal with heavy scalars, has 
been used. 
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Figure 16: The branching ratios for the radiative decay X2 ~^ Xil cis a function of n in the 
scenario 1 with M5 = 10^ GeV and 10^° GeV, and M2 chosen such that the rehc density is 
compatible with WMAP. The program SDECAY [17] has been used. 

virtuahty of the Z boson and by the additional Zff coupling. Furthermore, the Zx^x'j 
coupling needs to be strongly suppressed, thus the initial and final neutralinos need to be 
either pure gauginos or pure higgsinos. However, as discussed in section 13. 2[ for the LSP 
to form the DM in the Universe in such conditions its mass should be in the TeV range, 
see figs. [9]-[T0l Therefore, for relatively light neutralinos that lead to the cosmological relic 
density favored by WMAP (and, hence, have at least a small higgsino component leading to 
a non-negligible coupling to the Z boson), the branching ratio BR(x° — > Xil) is expected to 
be very small. 

This is exemplified in fig. [16] where the branching ratio for the decay of the next-to- 
lightest neutralino into the LSP and a photon, BR(x2 Xil)^ is displayed in the universal 
scenario 1 as a function of n for Ms = 10^ GeV and 10^° GeV; as usual, the value of M2 
is adjusted in such a way that the WMAP DM constraint is fulfilled. As can be seen, the 
branching ratio hardly reaches the level of 1% for Ms = 10^ GeV. For Ms = 10^^ GeV, the 
branching fraction is even smaller as the splitting M2 — Mi is larger, see table 1, leading to 
a more favored phase space for X2 ~^ Xi^* which controls the total decay width. 
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4.3 Gluino decays 

As already discussed in section 3.3, the gluinos decay through virtual heavy squark exchange 
either into neutralinos or charginos and a quark-ant iquark pair [46], g — > x'iQQ or g xfQQ'^ 
or into the two-body neutralino-gluon final state [49], g x^9i which is mediated by loops 
involving quarks and heavy squarks. Thus, the final state topology will consist of a neutralino 
or a chargino (which, if it is not the LSP neutralino, will subsequently decay according to 
the discussion held in the previous subsection) and one or two hard jets. 

It is experimentally important to know the number of final-state jets and, thus, the 
relative magnitude of the branching fractions for the-loop induced decays and the tree-level 
three-body decays. We have thus scanned the usual [M2, /i] parameter space and delineated 
the areas in which the branching ratio BR(^ — > g^iX^)i summed over all four neutralinos, 
is larger than 1,5, 10,25 and 50% (from lighter colors to darker ones). The results for the 
universal gaugino mass scenario are shown in fig. [T7| for a common scalar mass values of 
Ms = 10'^ GeV (left) and Ms = 10^° GeV (right); as usual the areas in which the WMAP 
constraint is satisfied and those excluded by the LEP bounds are also displayed. 

As can be seen, the branching ratio BR(^ g - Xi) is larger for lighter gluinos and thus 
smaller values of M2, and for larger values of the scalar quark mass M5. For instance, the 
g Xi branching ratio exceeds the level of 50% for M2 < 200 GeV and Ms = 10^° GeV. As 
discussed in Refs. [8,34], this is due to the fact that the main contribution to the radiative 
decay originates from loops involving top quarks and squarks (which have large couplings, 
oc rrit, to the higgsino components) and the ratio between this decay and the three-body 
decay scales as 

m'Jml [1 - log(M|K)]^ 

In particular, the branching ratio is substantial in the region where the cosmological relic 
density is generated by LSP annihilation through Higgs boson exchange, in which all neu- 
tralinos are relatively light and thus kinematically accessible. 

Similar results are obtained in the case where SUSY breaking occurs through an F-term 
that is not an SU(5) singlet. As an illustration, the areas of the [M2,/i] plane in which 
BR(^ — > g J2i Xi) is larger than 1, 5, 10, 25 and 50% are shown in the left-hand side of fig. [18] 
for the scenario 24 with Ms = lO'^ GeV. As can be seen, for the same value of M2 the 
branching ratio is significantly larger than in the scenario 1. For Ms = 10^^ GeV, which is 
not shown, the decay g ^ g Yli X^i is by far dominating compared to the three-body decay 
for low M2 values. The branching ratio for the string-inspired Oil model is shown in the 
right-hand side of fig. [181 In this case the radiative g ^ g Xi decay is dominant in a 
large band in which the mass of the gluino is close to that of the LSP neutralino, even for 
Ms = 10^ GeV. This band intersects the one in which the LSP has the correct density. 
Thus, there are sizable areas of the parameter space in which the radiative gluino decay is 
significant and even dominant. 
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Figure 17: The branching ratios for the decays of gluinos into neutrahnos and gluons, 
BR(^ —>■ gYliX^) the [M2,fj] plane in the scenario 1 with Ms = 10^ GeV (left) and 
10^° GeV (right); the region in which flh'^ is compatible with WMAP is also shown. 




Figure 18: The same as in fig. [T7] for the scenarios 24 (left) and Oil (right) with Ms 
10^ GeV. 
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5 Conclusions 



We have performed a comprehensive analysis of the MSSM in the scenario where all the 
scalars, except for the SM-like Higgs boson, are very heavy. This model, commonly known as 
Split Supersymmetry, is phenomenologically viable and much more predictive than the usual 
MSSM with light scalars. Indeed, besides the three soft-SUSY breaking mass parameters 
for the bino, wino and gluino, one has as basic inputs only the common scalar soft-SUSY 
breaking mass parameter M5, the higgsino mass parameter /i, which is not fixed by the 
requirement of proper radiative electroweak symmetry breaking, and the parameter tan/3. 
The model retains the interesting features of the usual MSSM, that is, it leads to a consistent 
unification of the gauge coupling constants at the GUT scale and provides a solution to the 
Dark Matter problem in the Universe. Nevertheless, as the scalars are extremely heavy, a 
large amount of fine tuning in the Higgs sector is present. 

In a first step, we presented the model and summarized our precise determination of the 
masses of the SM-like Higgs boson, the charginos, neutralinos and gluino, and of the different 
couplings of these particles. In our computation we resum large logarithmic corrections by 
means of the appropriate RG evolution, and we also include finite one-loop corrections to 
the particle masses. We performed an analysis of the variation of these masses with respect 
to the renormalization scale, which can be viewed as as rough estimate of the higher-order 
corrections. We then discussed the boundary conditions for the soft SUSY-breaking gaugino 
mass parameters and considered, besides the universal scenario in which they are set to a 
common value at the GUT scale, two representative sets of scenarios in which the GUT 
values are non-universal: one where SUSY breaking occurs via an F-term that is not an 
SU(5) singlet and another based on an orbifold string model. Finally, we briefly described 
how this model is implemented in the RGE Fortran code SuSpect which calculates the SUSY 
and Higgs particle spectra in the MSSM. 

In a second step, we summarized the available constraints on the model, first from collider 
searches and high-precision measurements and then from the requirement that the cosmolog- 
ical relic density of the lightest neutralino, which is expected to form the Dark Matter in our 
Universe, is compatible with the measurements made by the WMAP satellite. These analyses 
have been performed in the cases of universal and non-universal gaugino mass parameters. 
It turns out that several new features appear in the non-universal scenarios: in the case of 
the DM constraint, new channels such as neutralino annihilation through the exchange of a 
Z boson (which is ruled out by collider constraints in the case of universal gaugino masses) 
and gluino co-annihilation (which does not occur in the universal scenario) are possible. Fur- 
thermore, the annihilation of neutralinos through the exchange of a Higgs boson which then 
decays into a real and virtual W boson, leading to the three-body H —+ W j] final state, 
has to be taken into account. We also discussed the constraints on the common scalar mass 
parameter Ms from the requirement that the gluino lifetime does not exceed the age of the 
Universe, and showed that there are small differences in the non-universal cases compared 
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to the case of a universal gaugino mass parameter. 

We finally analyzed the phenomenology of the model, focusing on the decays of the Higgs 
boson and those of the charginos, neutralinos and gluinos. We emphasized the differences 
between the universal and non-universal cases and, for the Higgs boson decays, the differences 
between the SM and SUSY cases. For instance, we have shown that the invisible Higgs-boson 
decay H — XiXi can be substantial, reaching branching ratios of the order of 10% in models 
where the bino mass parameter Mi is much smaller that the wino mass parameter M2. The 
virtual contributions of charginos to the two-photon decay of the Higgs boson have also 
been shown to reach the level where they can be observable at the 77 option of the future 
linear e^e~ collider. Decays of heavier charginos and neutralinos into lighter ones and a 
Higgs boson have been shown to be potentially large, opening the possibility of measuring 
the Hxx couplings, while the loop-induced decays of heavier neutralinos into lighter ones 
and a photon are very rare, barely reaching the 1% level. Finally, we have compared the 
three-body decays of gluinos through the virtual exchange of heavy squarks, g — >■ qqx-, with 
the loop-induced two-body decays g gx, and have shown the latter to be potentially 
dominant in some non-universal scenarios. 
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Appendix 



In this appendix we present for completeness all the formulae for the one-loop radiative 
corrections that we include in the computation of the mass spectrum of the low-energy 
effective theory of Split SUSY. We largely follow the notation and the results of ref. [14], 
adapting when necessary the formulae of that paper to the case of Split SUSY. 



Gauge and Yukawa couplings 



We start by summarizing how the minimally renormalized gauge and third-family- Yukawa 
couplings, as well as the electroweak parameter v (from now on a hat denotes minimally 
renormalized parameters of the Split-SUSY effective theory), are extracted at the scale Q = 
Mz from the set of SM input parameters [as{Mz) , a{Mz) , Mz , Gp , Mf , M^- , mft(mb)]. The 
strong gauge coupling as is extracted from the SM input as{Mz) as 
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where 5s is a conversion factor that depends on the choice of the minimal renormalization 
scheme, i.e. 5^ = in MS and 5s = DR. Similarly, the electromagnetic coupling a is 

extracted from the SM input a{Mz) 
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where 5e = in MS and 5e = 1/3 in DR. The weak mixing angle (we denote sin^^y cind 
cosOw as s and c, respectively; we thus define = g''^ /{g^ + ^'^) and = 1 — M^/M"^ ) is 
computed from 
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where Uwwip"^) and Hzzip^) are the transverse (and real) part of the W and Z self-energy, 
respectively, and will be explicitly given below. The p parameter is defined as 
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while the quantity 5vB) which accounts for the vertex and box corrections to the muon decay, 
is: 
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In Split SUSY the non-SM contributions to 5vb are suppressed by the large sfermion masses 
and can be omitted. Once a and c^s^ have been computed by means of eqs. flA2p and flA3p . it 
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is trivial to extract the low-energy boundary conditions on the running electroweak coupling 
constants g{Mz) and g'{Mz). 

The top Yukawa coupling is extracted from the physical top mass Mt and the running 
electroweak parameter v according to: 

ht{Mz) v{Mz) = Mt + T.t{M^) , (A6) 

where Et(M|) is real part of the top quark self-energy (explicitly given below) computed 
with external momentum = M^, and -O is defined as: 

«-(Q)^2 ^| + "-<;f . (A7) 

In the low-energy effective theory of Split SUSY the couplings of the Higgs boson to the 
down-type fermions are SM-like, thus they do not have a large impact on the phenomenology. 
However, we include all the third-family Yukawa couplings in our analysis for completeness. 
The Yukawa couplings of the b quark and of the r lepton can be extracted from the running 
masses evaluated at Q = Mz according to: 

hAMz) v{Mz) = m,,r{Mz) + T^'^Mlr) , (A8) 

where T,f^^{ml^) contains the weak-scale threshold contributions to the b and r self-energies. 
In order to resum to all orders the potentially large logarithms of the ratio nib^r/Mz, the 
running masses frih^r are evolved up to Q = Mz by means of the RGE of the model with five 
quarks and SU{3) x t/(l)em gauge symmetry: 

m,,,{Mz) = mb,.(m,,.) 1 - ^ log ^ 1 _ log ^ , (A9) 

where {f3,, (3,) = (23/3,-80/9), (7,^7,^) = (8,2/3), and {Ys.Ye) = (0,6). The running 
bottom mass mb{m}j) is taken as input, while the running tau mass mT-{mr) is extracted 
from the physical mass Mr according to: 

m,(m,) =Mr (1- cr^ (AlO) 



where cr = A in MS and cr = 5 in DR. 



One-loop self-energies 

We now provide explicit formulae for the one- loop self-energies appearing in eqs. (]A3p . (]A4p 
and ( 1A6l) -( lA8l) . The transverse parts of the gauge boson self-energies Hvv (with V = Z, W) 
can be decomposed in a SM contribution Il^y and a chargino/neutralino contribution Hyy- 
In the Feynman gauge, that we adopt throughout this appendix, the SM contributions read: 
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16 vr^ 



W "^Z -2 



Bo(mw,mw) + 



+ {9u + 9%) H{mf, TUf) - Agj^ gj^ m) Bo{mf, m/)] , (All^ 

/ 



;i + 8c ) B22{mz,mw) 



8B22{mw^O)+^P^Bo{mw,0)] + ^ m^J .(A12) 



The Passarino-Veltman functions Bq, B22 and if are defined as in the appendix B of ref. [14]. 
The Higgs and gauge boson masses appearing in the equations above are interpreted as 
running masses. The summation in the last hne of each equation is over the fermion species, 
NI is the color number (3 for quarks and 1 for leptons) and gf = ll — ejs^ are the weak 
neutral-current couplings. and Avi^ are DR-MS conversion factors: they are both equal 
to zero in the DR scheme, while in the MS scheme A^ = — 2/3c^p^ and A^y = — 2/3j9^. 

The chargino and neutralino contributions to the gauge boson self-energies can be ex- 
pressed as: 

^, \ I ( 4- 1 Wfrr,^„,„.^^ 



16 7r2n|^(p^ 



Him^o, m^o) + 4 a°-^ 6°.^ m^o m^o Bo{m^o,m^o] 



+ 



J2[{''tjZ^ + KjZ^)H (m^f , ) + 4 6+^ Bo (m^+ , m^ 



(A13) 



16 7r2n^^(p2 



I (ajjw^ + hjw'^) H{m^o,m^+) + 4 ajj,y bijw m^o Eo(m^o, 

(A14) 



In general, we write the Feynman rule for the chargino or neutralino couplings to a gauge 
boson as —i 7^ [aPi + h Pr), the rule for the couplings to a scalar as —i {aPi + b Pr) and the 
rule for the coupling to a pseudoscalar as {aPi + h Pr)- Under this convention, the chargino 
and neutralino couplings to the gauge bosons are 



>'ijZ 



(A15) 



a i& - s^) 
^%z = 9cVaV,, + ^^ ' V,2V,2, 



btjz = 9cUaU,, + ^^^^^^U,2Uj2, (A16) 
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aijw = - 9cNi2Vji + -^Ni4Vj2 



-gcN,2U,i-^Ni^U,2 . (A17) 



The matrices and f/, V rotate the neutrahno and chargino states, respectively, so that the 
mass matrices N AiN N'^ and U A4cV^ are diagonal. The matrices AiN and Aic are given 
in eq. 1^. We assume that there are no CP- violating phases in the higgsino and gaugino 
mass parameters, and we choose A^, U and V to be real, allowing for negative signs in the 
chargino and neutralino masses. 

The top quark self-energy in eq. ( ]A6I) can be expressed as 



St mj 

nit 



\ 4e 



3 log - c/j ) + 



9 



cr 



+ 



4 + hi 



Bi{mt,mH) + Bo{mt,mH) +Bi{mt,mz) - Bo{mt,mz) 



Bi(mb,mw) 



h 



9' 



J Bo{mb,mw) + 2 \_ 



9' 



Bi(mt,mz) + S_ 



R 



4^9tL 9tf 



Bi{mb,mw) + 5r 
Bo(mt,mz) + Sr 



(A18) 



where {cr,6r) are equal to (4,-1/2) in MS and to (5,0) in DR. The Passarino-Veltman 
functions Bq and Bi are defined as in the appendix B of ref . [14] . The weak-scale contribution 
to the bottom quark self-energy T,f^{mb), appearing in eq. (lASp . can be extracted from 



eq. ( ]A18I) by omitting the first line and replacing everywhere t ^ h. The analogous quantity 
for the tau lepton, S*^(m^), can also be extracted from eq. ( 1A18I) : one has to omit the first 
line and the terms proportional to hi in the third line, then replace everywhere t t and 
TRb 0. 



Corrections to the Higgs mass 

We provide here the formulae for the Higgs mass corrections 5^^{Q) and 5'^{Q) appearing 
in eq. (JTj). The SM contribution reads [24] 

Gf Ml 



5^^(g) 



e/i(e,Q) + /o(e,Q) + rv-i(e,Q) 



v/2 16 7r2 

where ^ = mjj/Ml and the functions fk{^, Q) are defined as: 



(A19) 



/o(e,g) 



(A20) 



-6 log 
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^iog4 + 2z(r^) + 4c^z(c^r^) 
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H 



4 log 
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Z L 



:i + 2cV24^ 



Mi 



log 



6Z(r') - 12 c^Z(c2r') - 12 logc^ 



mi 
Ml 



2 + Z 



7717 



m 



(A21) 



(A22) 



Z(x) 



In the equations above we define = M^^/M^ in terms of physical masses. On the other 
hand, as described in section 2.2, we explore the consequences of choosing as either the 
physical or the running top mass. The auxiliary function Z appearing in the equations above 
is defined as: 

2Atan-i(A-i) (^>i) 

(A23) 

Alog[{l + A)/{l-A)] {x<l) , 

|l-4x|V2. 

The results of ref. [24] are derived under the assumption that the Higgs quartic coupling 
A appearing in eq. ([7]) is expressed in the MS renormalization scheme. If A is expressed in 
the DR scheme the Higgs mass correction in eq. ( 1A19P is modified as 

^2 /if2 



where A 



16 vr^ rn^j 



2c' 



(A24) 



The chargino and neutralino contributions to the Higgs boson self-energy and tadpole, 
appearing in the correction term 5^{Q) defined in eq. (IE]), read 



(A25) 



(A26) 

The Passarino-Veltman functions G, Bq and Aq in the equations above are defined as in the 
appendix B of ref. [14]. The chargino and neutralino couplings to the Higgs boson are 



(A27) 
(A28) 



where gu,d and g'^ ^ are the effective Higgs-higgsino-gaugino couplings defined in eq. (P). We 
define S^ki^ij} = SkiSij + 6kjSii, where 6ij is the Kronecker delta and summation over repeated 
indices is understood. 
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Corrections to the chargino and neutralino masses 

We provide liere tlie formulae for tlie radiative corrections to the chargino and neutralino 
mass matrices, once again adapting to the Spht-SUSY case the results of ref. [14]. The 
one-loop neutralino mass matrix reads 



where 



The one-loop chargino mass matrix is instead 



Mcip') =M'c- Kip') Mc - Mc ^tip') - ^sip') 



(A29) 
(A30) 

(A31) 



The tree- level mass matrices M% and A4q are given in eq. and are expressed in terms 
of minimally renormalized parameters. The neutralino and chargino self-energies Hi glp"^) 
and S^^(p^) are 4x4 and 2x2 matrices, respectively, and they read 

2 

1 — 1 



k=l 

4 



16 7r^S°., 



^ 4 
k=l 

2 

= '^J2^xt{ «ifcG+ Bo ,mw)-4: a%w 



(A32) 
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Bo{m +,mw 



+ Sr] } 



k=l 



^ "ikZ "'jkZ 
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(A33) 



16 7r^S+^. 
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HiG+ KjG+ Bi{m^o , mw) + 2 a^.,!^ a 



''kiW "'kjW 



Bi{m^o.,mw) + Sr 



k=l 
2 



fc=i 



Bi{m^+,mz) + Sr +2a^^a+fc^ Si(m^+, 0) + 5/? 



+ 2 o-ifcz ^jfcz 



(A34) 
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fc = l 

2 



kiW ^kjW 



m 



k=l 



4 ^ifcZ '^jTcZ 



Bo{m^+,mz) + Sr 



4ft+ fl+ 



(A35) 



where 6r = in DR and 6r = —1/2 in MS. The formulae for the self-energies S^j and 
can be obtained from those for S° and S^, respectively, by replacing a — > 6 in eqs. flA32p 
and (lA34p . In the equations above we denote by ai-,$ and byq, the couplings of a bosonic 
field $ = (Z, VF, 7, H, (7°, G+) with one rotated neutralino (or chargino) mass eigenstate and 
one unrotated neutralino (or chargino) gauge eigenstate. 

In particular, the couplings of an unrotated neutralino ip^, a rotated chargino and a 
charged pseudo-Goldstone boson read 



-la. 



jk 



Vjk { SiiSk2 + 5i25k2 + ^i44l 



(A36) 
(A37) 



The couplings of an unrotated neutralino ip^, a rotated chargino x'j ^ W boson read 

1 



g Vjk ( -Si2Ski + -j= Si4,6k2 



jk 



(A38) 
(A39) 



The couplings of an unrotated neutralino ipi, a rotated neutralino x'j aiid a neutral pseudo- 
Goldstone boson read 



«?;G0 = - b^ijGO = 2 [~9'd ^{il^kS} - g'u ^{il^ki} + gd S{i2Sk3} + (ju 5{i244} j • (A40) 

The couplings of an unrotated neutralino ip^, a rotated neutralino and a Higgs boson read 



HjH 



%H = ^ g'd S{iiSk3} + g'u ^{ii44} + gd S{i25k3} - gu ^{^244}) • (A4l) 



The couplings of an unrotated neutralino ip^, a rotated neutralino Xj and a Z boson read 

g 



2c 



Nik. 



(A42) 
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The couplings of a rotated neutralino x^, an unrotated chargino ipj' and a charged pseudo- 
Goldstone boson read 

dt^+ = Nik (^-^ SkiSj2 + (^k2Sj2 + Qu Ski^ji^ , (A43) 

K^G+ = (j^SkiS,2 + ^Sk2S,2-~gdSk3Sjiy (A44) 
The couphngs of a rotated neutrahno Xi", an unrotated chargino ipj' and a W boson read 

= 9 Nik (^-Sk2Sji + Sk^dj^ , (A45) 

bfjw = 9 Nik (^-Sk2Sji - -^Sk3Sj2y (A46) 

The couphngs of an unrotated chargino ipt^ ^ rotated chargino and a neutral pseudo- 
Goldstone boson read 

a^jGO = (gd 5iih2 - (Ju Si2Ski^ , (A47) 

b^jGO = ^Jk (qu kl5k2 - Qd <^j24l) • (A48) 

The couplings of an unrotated chargino ipt^ ^ rotated chargino Xj' and a Higgs boson read 

afjH = -^Ujk(^gdSiiSk2 + guSi2Skij, (A49) 

KjH = -^^jfe(^«^il5fc2 +^d<5i24l)- (A50) 

The couplings of an unrotated chargino ip^y ^ rotated chargino Xj' and a Z boson read 

SiiSki + ^ ^i242 j , (A51) 

^2 "2 

fcj^ = g c Ujk (Sii6ki + ^ ^^242) ■ (A52) 

Finally, the couplings of an unrotated chargino tpf, a rotated chargino x'j and a photon read 

e Vjk (Sii6ki + 5i242) , (A53) 



= eUjk(^SiiSki + Si2Sk2y (A54) 

The self-energies in eqs. ( ]A29p - (]A3ip induce an external momentum dependence in the 
neutralino and chargino mass matrices. We compute each of the physical neutralino and 
chargino masses by diagonalizing the one-loop mass matrix with the external momenta in 
the Passarino-Veltman functions set equal to the tree-level mass of the corresponding particle. 
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More one-loop translations from DR to MS 



Beyond tree level the boundary conditions on the quartic Higgs coupling and on the Higgs- 
higgsino-gaugino couplings given in eqs. ©-([S]) are valid only in the DR renormalization 
scheme. The MS scheme breaks supersymmetry, therefore in that scheme the SUSY rela- 
tions between the gauge couplings and the gaugino couplings and those between the gauge 
couplings and the quartic scalar couplings are not preserved [50]. In the MS scheme the 
boundary conditions in eqs. ([3])-([S]) are modified as 
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(A55) 
(A56) 
(A57) 
(A58) 
(A59) 



where AthA is the 0{h^) threshold correction defined in eq. (E]). 

For completeness we conclude this section by providing the relations between the DR 
and MS definitions of the gaugino and higgsino mass terms: 



MS 



DR 



(A60) 
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